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Abstract. Universal properties of efficiency at maximum power are in¬ 
vestigated in a general setting. In particular, it is demonstrated how 
successive symmetries placed upon the dynamics manifest themselves 
at the macroscopic level. A general condition is derived for which ther¬ 
modynamic devices are able to attain a reversible operation. 


1 Introduction 


According to Callen [T], thermodynamic theory is rooted in the statistical properties 
of large systems and the symmetries of the underlying fundamental laws of physics. 
A profound illustration of the latter was given in 1931 by the seminal work of Lars 
Onsager [5], demonstrating that time reversibility of the microscopic dynamics be¬ 
comes evident at the macroscopic level through reciprocity relations among the lin¬ 
ear response coefficients. Experimental evidence of such relations was already known 
to exist in various thermoelectric systems, cf. the Seebeck and Peltier effects. With 
the development of stochastic thermodynamics |3|4|5 |6|7|8j , such reciprocity relations 
were shown to exist also among the nonlinear response coefficients, see for example 
pru] . These relations are a consequence of the fluctuation theorem, entailing a sym¬ 
metry property for the probability distribution of entropy production. 

Entropy production plays an important role in the performance of thermal engines. 
Stochastic thermodynamics provides the natural language to describe such systems. 
Recent work has concentrated on the efficiency of small sized thermodynamic ma¬ 
chines, with a focus on the efficiency at maximum power (EMP) output 
[1 1 n 211,311411511 fill 7118119120121 l22l23tMl^ . Historically, EMP for heat engines was 
first discussed by Novikov [55] who obtained the following result: 


f]=l - ^J\ - T]c = 


)7c 

2 8 



( 1 ) 


with rjc the Carnot efficiency. This result was later rediscovered by Curzon and 
Ahlborn m- Quite remarkable, the resulting EMP is device independent and con¬ 
tains only the temperature ratio. This raises the question whether this efficiency can 
be placed at the same level as the Carnot efficiency. As was shown in m, in the 
case of heat engines for which the heat and work flows are proportional, the first 
order term of Eq. [^is recovered. When an additional symmetry is present, also the 
second order term was reproduced in [T^ . However, calculations of the EMP for var¬ 
ious model systems have revealed a large variety of EMP expressions. The purpose 
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Fig. 1. The generic setup consisting of an arbitrary mesoscopic system in simultaneous 
contact with different reservoirs. Each reservoir is characterised by the values of the intensive 
parameters such as temperature, pressure or chemical potential. 


of this work is to investigate EMP in a model independent setting and for general 
energy conversion machines. We demonstrate how symmetries and constraints at the 
microscopic level, combined with the fluctuation theorem, emerge at the macroscopic 
level via the expression for the EMP. 


2 A generic setup for energy conversion 

Our starting point is the generic setup of a system in simultaneous contact with 
different reservoirs as shown in Fig. Each reservoir is characterised by the values 
of its intensive variables such as temperature, pressure and chemical potential. When 
these values are different for the various reservoirs, irreversible processes occur trying 
to establish a common final equilibrium state. These processes involve the exchange 
of extensive quantities like energy and particles for example. We consider infinitely 
large reservoirs so that a non equilibrium steady state is reached. For concreteness, 
we consider a setup which allows for the exchange of two extensive variables Xi and 
X 2 when the associated thermodynamic forces (affinities) Ai and A 2 are nonzero. 
These forces are defined as the difference of intensive variables, 

Ak=Fl,-Fk, ( 2 ) 

which in turn are determined from the entropy fundamental relation S'dXi}) as Fk = 
dS/dXk n. Such a setup with two forces is natural in the context of thermodynamic 
engines where one force (e.g. a temperature difference) provides the driving energy 
to do work against the other force (e.g. a chemical potential difference). 

The exchange processes are strongly determined by the details of the intermediate 
system and by its connection to the reservoirs. They are however not completely 
arbitrary, but must be in agreement with the microreversibility of the underlying 
microscopic dynamics. Such properties of the exchange processes can be quantified 
through the use of the cumulant generating function (CGF). Denoting Pt{Xi,X 2 ) as 
the probability to have a net exchange of Xi and X 2 for the two extensive quantities 
during the time interval from 0 to t, the CGF is defined as 

1 1 A ’, 

G{X-A) = lim - ln(e-("^^^+"=^^)) = lim - ^ ^—^{{XlXi)}X\Xi. (3) 

>-oo t t^oo t z!?! 

i,j=0 

The cumulants ((XJX^)) represent measurable macroscopic quantities and are deter¬ 
mined by the specific system and its connection to the reservoirs. In particular, they 
depend on the intensive parameters. Replacing by F^. +Ak and expanding in terms 
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of the forces gives: 

OO 

\^{{XlXi))/t= Y, ( 4 ) 

k,l—0 

The response coefficients are determined with respect to the equilibrium state 
corresponding to F'f. = F^- Time reversibility imposes the following symmetry relation 
upon the generating function 

G{X;A) = G{A- X]A). (5) 

As a consequence, after substituti ng p| into one can make, by rearranging the 
terms, the identification (see also Ml: 


^ ^ f 1 \i-\-j-\-m-\-n 

jkl _ V j-k—m,l — n 

m—0 n—0 


( 6 ) 


This result is central to the discussion below, and expresses the macroscopic mani¬ 
festation of time-reversibility at the microscopic level. 

For energy conversion, the focus is on the first two cumulants, i.e. the time averaged 
currents 

Ji = lim {{Xi))/t and J 2 = lim {{X 2 ))/t. (7) 

>-oo >-oo 

During steady state operation these currents, together with their corresponding forces, 
determine the entropy production rate via the bilinear form [1]: 


S — J\A\ -j- J 2 A 2 ^ 0. 


( 8 ) 


In the context of thermodynamic engines, this relation can be used to obtain a general 
and device independent definition of efficiency. Without loss of generality we consider 
J 2 A 2 > 0 so that the flow J 2 is spontaneous and regarded here as the driving process. 
A 2 is the corresponding driving force. This process is used to induce a flow Ji against 
the force Ai, so that JiAi < 0. The efficiency of the conversion process can be defined 
as follows: 


0 < 77s 


—AiJi 

A2J2 


< 1 . 


(9) 


This concept was coined entropic function in |31j . 


3 Microscopic constraints 

From Eq. it is straightforward to show that Onsager symmetry = L}^ is 
present at tne linear order in the currents. Combining Eq. with Eq. Q leads to: 

Ji = (L°o/2) Ai + LA 2 + (T 20 / 2 ) Al+[Q + L°l/2] A 1 A 2 + MAl + ... (10) 
J 2 = LAi + {Lll/2) A 2 + QAl + [Lll/2 + M] A 1 A 2 + {L°1/2) AI + ... (11) 

with 

L = L0?/2 ; Q = L)?/2 ; M = 1^1/2. 

More surprisingly, a similar structure is also appearing at higher orders. 


(12) 
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3.1 Strongly coupled flows 


In their paper EH, Kedem and Caplan introduced the concept ’’degree of coupling” 
between the flows Ji and J 2 . Since their work was limited to the linear order terms 
in Eqs. (10) and (11), the degree of coupling q was defined in terms of the linear 
coefficients q = L/ a/LJqLq}. Maximal coupling is achieved for g = ±1, in which case 
the two flows are proportional to each other. 

The concept of maximal coupling can be extended beyond the linear regime as follows. 
We call a process strongly (maximally) coupled if the exchange between reservoirs is 
constrained by the property: 

= eX2. (13) 


This property will be further referred to as strong coupling of the process (SCP), and 
is especially relevant in nanoscale systems and models [T6l3^ . Obviously it implies 
Ji = eJ 2 - At the level of the response coefhcients it entails: 




k,l 


Combining this with the result from the fluctuation symmetry gives: 


(14) 


Ji — L + A.2) + (eA-i + A2) (QAi + MJ\. 2 ) 

+ (eA-i + . 4 . 2 ) (i?A^ + S'w4.iw4.2 + ... (15) 

One notices immediately the appearance of the combination sAi + A 2 . A straight¬ 
forward calculation, given in appendix shows that this is true at every order. 
Hence, one can identify a unique stalling force Ai = —A^jE for which both fluxes 
simultaneously vanish. The efficiency is then 


—AiJi — eAi 

A 2 J 2 A 2 


(16) 


Hence SCP is a sufficient condition to attain reversible operation at the stalling force. 
Notice also that this result is independent of the specific value of e. 


3.2 Strongly coupled forces 


In |15j a stronger condition was considered. When the generalised forces only appear 
in the combination eAi + A 2 , the probability distribution Pt,yti,^ 2 (Ai, A 2 ) can be 
written as Pt,eAi+A 2 {^ii^ 2 )- We will refer to this property as strong coupling in the 
forces (SCF). As consequence, the generating function has the property 


dG{X-A) JG{X;A) 
dAi ~ ^ dA2 


which can be used on the definition of the generating function Eq. ^ to derive a 
relation similar to Eq. (14) 


(fc + i)4f’' = e4j+i(; + i). (18) 

Using this relation, it is easy to see that we have a collapse into a single generalized 
force A = eAi + A 2 , thus also making the current expansions only dependent on A, 
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with only one free coefficient at any order. Indeed, successively applying Eq. (18) on 
yields: 




Hence 


(19) 

( 20 ) 


>^1 = E E = E + a ^ y . 

p=0 i=0 p=0 

For the determination of the second current J 2 , we show in Appendix B that the 
combination of SCF and the fluctuation symmetry automatically implies SCP, that 
is Ji = eJ 2 - The proportionality constant between the currents e is the same as in 
eAi + A 2 ' 

Ji=LA + MA^ + NA^ + 0{AY-, (21) 

J 2 = Ji/e. (22) 


4 Implications on the efficiency at maximal power 

It is clear from the previous discussion that the number of kinetic coefficients is 
reduced when microscopic constraints are imposed upon the system. We now clarify 
how this reduction of coefficients is reflected in the efficiency at maximum power. The 
maximum is determined in the following way: given the driving force A 2 , what value 
of the loading force Ai maximises the power output —AiJi? The answer requires 
solving the following nonlinear equation: 

5 . . _ . . dJi 


dAi 


(— AiJi) — —Ji — A 


^dAi 


= 0 . 


(23) 


Other optimisation schemes (see for example EEI]) are not considered here. 

Since an analytical solution is not possible in general, we resort to a series solution 

Ai = C 1 A 2 + C 2 .A 2 + C 3 A 2 + ... (24) 

and d etermi ne the coefficients. We already gave the current in its most general form 
Eqs. ( lOpT ), due to the FT. Solving for maximal output and using the coupling 
parameter q, the EMP to first order in A 2 reads 


q 


V = 


v/ZM (9^(4 + q^)LW2 + (8M + q^Ll°)) 

- 4 LOO 9 ((4 + g2) Q + 2L%) L™ - 8g(L™)2L0i 


’A 2 


4-2g2 ’ 8iL^YL^iq^-2Y ^.,...(25) 

For a system with SCP, the coefficients from Eq. ( [24| appear directly in the efficiency: 

■q = -£AiIA 2 = -£ (ci + C 2 A 2 + C 3 A 2 + ...). (26) 

The efficiency expansion substantially simplifies: 

1 QA 2 

JLe 


V = 


_ 2LR - 2MQe + 2LS£ 2 


16L2£2 


(27) 


From which we already see at lowest order the known universal factor 4 HU. The 
latter expansion is especially relevant for autonomous isothermal motors since the 
same definition of efficiency m applies. The expansion simplifies further when the 
forces are strongly coupled (18). From Eq. (27) and substituting Q = eM, R = e^N 
and S = 2eN we find: 




1 M , 

— — — — A2 

2 8L 


m2 - 2LN 




OiAl 


16L2 


(28) 
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5 Discussion and perspectives 


The transition of the 18th to the 19th century marks the beginning of the industrial 
revolution. A prominent ingredient was the development of efficient steam engines. 
James Watt and his colleague engineers succeeded in improving the performance 
by a factor of twenty. In the wake of these developments, Sadi Carnot published a 
manuscript entitled ” Sur la puissance motrice du feu”. He showed that the efficiency 
of thermal machines is bounded by that of a machine operating reversibly. This in¬ 
sight will lead, with the notable input of Clayperon, Clausius and Kelvin, to the 
development of thermodynamics. Clayperon drew the attention to the fundamental 
significance of the Carnot cycle. Clausius introduced the state function entropy and 
formulated the second law of thermodynamics. The name of Kelvin is attached to the 
scientifically well defined concept of temperature. 


In its original formulation, thermodynamics describes macroscopic systems at 
equilibrium. With the work of Onsager and Prigogine, it was extended to describe 
systems in local equilibrium. More recently, one has been able to formulate thermody¬ 
namics for small system far from equilibrium. Stochastic thermodynamics is arguably 
the simplest such formalism. The most spectacular advance are probably the integral 
and detailed fluctuation theorems, which replace the positivity of the entropy pro¬ 
duction by a symmetry property for its probability distribution. This result is valid 
for any size of the system and for any nonequilibrium state. Revisiting the question 
of efficiency of thermal and other machines from this point of view has led to some 
remarkable discoveries. Onsager symmetry implies that efficiency at maximum power 
with respect to the load in the regime of linear irreversible thermodynamics is at 
most half of the reversible efficiency HH, cf. Eq. ( |M| . The constraint, imposed by the 
fluctuation theorem, has also an impact on the nonlinear response coefficients |30j . 
When asking the right question about efficiency, it leads to universal values for the 
quadratic nonlinearity for strong coupling and an appropriate additional property. In 
the present paper, we have shown that strong coupling with fluctuation symmetry 
implies an expansion of the form given in Eq. (21). The additional requirement to fix 
the value of M depends both on the system under consideration and on the definition 
of the efficiency. For example, for appropriate systems with ’’left/right” symmetry, 
the coefficient M becomes zero [12]. For appropriate thermal machines with the effi¬ 
ciency ’’defined a la Carnot”, the coefficient is equal to 1/8 |TS|. These universal values 
can be seen as manifestations of micro-reversibility (Liouville’s theorem or unitary 
evolution). 


We finally mention another extremely exciting recent development. The efficiency 
of small scale systems is fluctuating. The fluctuation theorem has implications on the 
probability distribution for this so-called stochastic efficiency [33134135136157155] . The 
most surprising result is that the reversible efficiency becomes the least probable in the 
long time limit for time-symmetric driving |33l37j . For time-asymmetric driving the 
probability distributions (or more precisely the large deviation functions) for forward 
and backward driving cross at the reversible efficiency [34135] . Hence the reversible 
efficiency is a special point for stochastic efficiency. It remains to be seen whether 
this result has practical implications, with in particular the absolute measurement 
of the temperature. The situation is reminiscent of the Jarzynski equality, which 
allows the measurement of equilibrium free energy differences by a large number of 
nonequilibrium experiments instead of a single reversible one. Here the reversible 
measurement of Carnot efficiency is replaced by multiple irreversible measurements 
from which the least likely efficiency can be identified. 
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A Unique stalling force 

Using the fluctuation relation Eq. ([^ and the SCP Eq. (141, we will now give the 
proof of Eq.(15). The first step is to show that the following quantity 


A{p)=^{-eyLl:^^^ = Q. 


(29) 


z=0 


Applying the fluctuation theorem Eq. ml on the coefficients yields: 


A{p) = J2{-sy 


2 = 0 


p—t I 


_ _ ( I'll+m+n 


_m—0 n=0 


min! 


—2—m,2—n 
m+l,n 


(30) 


Repeatedly using SCP Eq. ( |14[ ) gives after n times: 

p 


A{p) = J2{-sy 
2 = 0 
P 

= 


l'_i'il+m+n 1 

V j p—i—m,i—n ^ 

m+l+n.O ^ 


2=0 


EE 

m=0 n=0 
2 p —i-\- n 

E J- X '' V Y p—i—{m —n), 2 —n 

(m' — 


n\e 

n=0 ra'—n 


(TO'-n)!^""'+i’° 


(31) 


In the second step the substitution m' = m + n was used. Next we substitute the 
index n = 0...i hy I = i — n = 


A{p) = Y,{-ey 


i=0 


1 (- 
(i — (m! 


(- 1 ) 






— l—m' ,l 


^ (rn'-(i-l)y. ™'+b0 

L/=0 ^ ' m'=i-l ' ^ '' 


(32) 


Now we simplify with the following substitution m" = m! + 1. The double primes are 
dropped for ease of notation. 


A{p) = Y^{-ey 


2 = 0 


2 ^ P 


(- 1 ) 


1+m—Z 


Ll 


— m,l 


,Z=0 


(i-iy.e^-i^ (m-iy. 

N / fm — 1 '■ ' 


(33) 


Next we take the coefficient j outside the summation over i. This requires a double 
switch of the index i. First the sums over i = 0...p and m = i...p can be replaced by 
the sums m = 0...p and i = 0...m where the order of summation has changed: 


p m I 


= E E E 

m=0 2=0 Z=0 


1+m —Z 


(-i)v (-1) 

{i — iy (m —i)! 


—m,l 


_ tP 

_ qi 


(34) 


Second, a similar switching is done from i = 0...m and I = 0...i to * = l...m and 
I = 0...m which makes 




1+m—Z 


m=0 Z=0 


(- 1 )* 


(i — l)\(m — iy. 


(35) 


The sum appearing between square brackets yields (—l)*dz,m and hence 


A{p) = = -A{p) = 0. 


(36) 


m=0 
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It can be used to find the general stalling force of the process. Using Eq. (29) on the 
first coefficient g of the p-th order in the current expansion of Ji gives 


jw = ^ - { sta )- ( 37 ) 

i —0 z—1 

Finally, since x* — y* = (x — 2/)(X]}=o this result proofs that every order p 

contains the factor sAi + A 2 - 


B SCF implies SCP 

We now show that strong coupling in the forces implies also strong coupling in the 
process. We will proof that when the thermodynamic forces Ai are strongly coupled 
Eq. (18), it automatically implies strong coupling in the process Eq. (14). In the main 
text it was already established that the strong coupling relation Eq. (18) allows us to 
write the expansion of the cumulant i, j as follows: 


{{XlXi))=tY,L°;fA^ 

p—0 

with A = eAi + A 2 - So to prove the strong coupling in the process, the relation 

tO,P _ ^tO,p 


must hold. From the strong coupling in the forces Eq. (181 it follows: 

Using the fluctuation relation Eq. (§ on the Ihs leads to 

rl,P_^^ (_l)b+^-+-+") - 

~ 2-^ 2-^ ■ ■ 


71^ 


m—0 n—0 


= E 


(^ — ly+j+n 


n—0 


ni 


- l,p—n 
''z,j'+n 






E 

n—0 


n\ 


- 0,p—n 
^2+1,j+n 


The last term can be rewritten with the fluctuation theorem Eq. (§: 


rO.p _ 


P ^_2^^i+l+j+n 


7.1 


r 0,p — n 
■^2+1,j+n' 


n—0 


So Eq. (41) can be written as 


= E 


(-l)^+^'+%l.p-n , rO,P 


n—0 


n: 




Rearranging the terms then leads to 


rO.p _ 


(p+i)L":r‘ - E +Dd-jv-r- 

n—0 


(38) 


(39) 


(40) 


(41) 


(42) 


(43) 


(44) 
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where Eq. (181 was used on the coefficients " 
then on 


of the last term. Using the FT Eq. 


tO,p 


= e 


p+i 


(P + 1) 

n=0 


(^ — ly+j+n 
nl 


j 0,p+l —n 


The two sums can be merged together: 


E 


n—O 


(^—\y+j+n 

n\ 


{p-n + l)L 


0,p—n+1 

2,j+n 


(45) 


rO.p 


P+1 

= ^E 


n=0 


n! 


{n)L 


0,p—n+1 


(46) 


The term n = 0 can be dropped and substituting the index m = n — 1 yields the 
desired result. 


rO,p 

^i+i,i 


= £ 


E 

_m=0 


(- 1 ) 


2+j+m+l 


ml 


j 0,p—m 
j+m+1 


= sL 


0,p 

*,i+i 


(47) 


where in the last step the FT Eq. § was used again. So the strong coupling in the 
forces has been shown to imply the strong coupling of the process Eq. (14|. 
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